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Abstract 

The Standard Model Higgs boson is searched for in 233.2 pb -1 of data collected 
by the L3 detector at centre of mass energies from 192 GeV to 202 GeV. These data 
are consistent with the expectations of Standard Model processes and no evidence 
of a Higgs signal is observed. A lower limit on the mass of the Standard Model 
Higgs boson of 107.0 GeV is set at the 95% confidence level. 
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1 Introduction 



The Standard Model of the electroweak interactions |T] contains a single Higgs doublet |2j which 
gives rise to a neutral scalar particle, the Higgs boson. Its mass, mn, is a free parameter of 
the theory. A global fit to electroweak precision measurements results in an upper limit on 
mR- of 133 GeV [Sj at 95% confidence level (CL). Results of SM Higgs boson searches in e + e~ 
collisions were published up to centre of mass energies of 189 GeV in the Higgs mass range up 
to 95.3 GeV by L3 J3] and by other LEP experiments [3]. 

In this paper, the results of a Higgs search performed on the data sample collected by L3 
at y/s up to 202 GeV are reported. The dominant Higgs production mode is 

e+e-^Z*^HZ . 

The processes of W + W~ and ZZ fusion, which contribute with smaller rate to the Higgs pro- 
duction in the Huu and He + e~ channels, respectively, are also considered. All significant Higgs 
decay modes are considered in the search. The largest sources of background are four-fermion 
final states from W and Z pair production, as well as e + e~ — >qq(7). 

2 Data and Monte Carlo samples 

The data were collected using the L3 detector |Hj at LEP during the year 1999. The integrated 
luminosities are 29.7 pb" 1 at ^ = 191.6 GeV, J 83.7 pb" 1 at y/s = 195.5 GeV, 82.8 pb" 1 at 
y/s = 199.5 GeV and 37.0 pb" 1 at = 201.8 GeV. 

The Higgs production cross sections and branching ratios are calculated using the HZHA 
generator [Jj. Efficiencies are determined on Monte Carlo samples of Higgs events, generated 
using PYTHIA [8J. Standard Model background estimates rely on the following Monte Carlo 
programs: PYTHIA (e+e" qq( 7 ) and e+e" -> Ze+e"), KORALW (e+e" W+W"), 
KORALZ PS] (e+e-^r+r-), PHOJET [TT] (e+e" ^e + e"qq) and EXCALIBUR [H] for other 
four fermion final states. The number of simulated events for the dominant backgrounds is at 
least 100 times the number of collected data events for such processes. For the Higgs signals, 
at least 2000 events are simulated for each search channel and for several masses. Higgs events 
are simulated with mn between 95 and 110 GeV, with a step of 1 GeV. Events for Higgs masses 
between 50 and 95 GeV are simulated with a step of 5 GeV. 

The response of the L3 detector is simulated using the GEANT program pQI], taking into 
account the effects of multiple scattering, energy loss and showering in the detector. Hadronic 
interactions in the detector are modelled using the GHEISHA program ^3]. Time dependent 
detector inefficiencies, as monitored during the data taking period, are also simulated. 

3 Analysis procedures 

The search for the Standard Model Higgs boson is based on the study of four distinct event 
topologies representing approximately 98% of the HZ decay modes: qqqq, qqw, qq£ + £~ (£ = 
e, /i, r) and r + T~qq. With the exception of the HZ^r + T~qq decay mode, all the analyses are 
optimised for H — > bb decay. This mode represents about 82% of the Higgs branching fraction 
in the mass range of interest. 

All these channels are analysed in three stages. First, a high multiplicity hadronic event 
selection, preserving most of the Higgs signal, is applied to reduce the large background from 
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two-photon processes. In a second stage, a cut based analysis is applied to the qqqq and lepton 
topologies, while a neural network based analysis is used for the qqz^ final states. All the 
analyses use topological and kinematical variables, which are not strongly dependent on the 
Higgs mass. B hadrons are identified on the basis of an event b-tag variable, obtained as a 
combination of the b-tag for each hadronic jet A neural network is used to calculate the 
b-tag for each hadronic jet from the three-dimensional decay lengths, properties of semileptonic 
b decays and jet-shape variables. The third part of the analysis is the construction of a final 
discriminant for each topology. This is built from a combination of the event b-tag variable and 
the reconstructed Higgs mass, for the cut based analyses. For the neural network based analysis, 
it is a combination of the neural network output with the reconstructed Higgs mass. The 
distributions of the final discriminants are computed for the data, for the expected background 
and signals for each Higgs mass hypothesis. 



The selection aims to single out events with four jets, two of which contain b hadrons, while 
the invariant mass of the other two must be consistent with the Z mass, mz- Background from 
Standard Model processes comes mainly from qq final states with hard gluons, W + W~ and ZZ 
events, especially those where one of the Z bosons decays into b quarks. 

A preselection to accept high multiplicity hadronic events is applied. The retained events 
are forced into four jets with the DURHAM algorithm ^B]. Then a kinematic fit requiring 
4-momentum conservation is performed. 

In order to differentiate the Higgs signal from background, selection criteria mainly based on 
kinematic variables, are chosen to maximise the expected performance of the analysis [3] . The 
quantities used for the selection are: the dijet masses, the minimum jet energy, the maximum 
energy difference between any two jets, the parameter of the DURHAM scheme for which the 
event is resolved from three jets into four jets, Y34, and the number of charged tracks. The 
consistency of the dijet masses with a given mn hypothesis is quantified by a x 2 -probability 
that depends on wh and mz [31 ■ Only a loose cut is placed on this variable. 

After these cuts, about 85% of the expected background comes from W + W~ events. These 
are characterised by their low b-tag values and by dijet masses close to the W mass. The 
remaining events are then split for the final analysis into a high purity and a low purity sample 
according to the value of the reconstructed dijet mass, M^. This is calculated with a five 
constraints kinematic fit where equal masses are assumed for the two dijet systems. The event 
is assigned to the high purity sample if M^P > a x m^ + b, or to the low purity sample otherwise. 
The quantities a and b are optimised for each centre of mass energy. Typical values are a ~ 0.6 
and b ~ 40 GeV. The low purity sample contains most of the properly reconstructed W + W~ 
events. Loose b-tag cuts are applied in the final selection and they are separately optimised at 
each centre of mass energy for the two samples. They are tighter for the low purity one. 

The low purity sample contains 109 events with 108 expected from background, correspond- 
ing to an efficiency of about 28% for HZ^bbqq and =105 GeV. The high purity sample 
contains 32 events with 40 expected from background, corresponding to an efficiency of about 
30% for a Higgs mass of 105 GeV. 

The final discriminant is then calculated as the weighted probability [3] that an event is 
consistent with the background distributions of both the b-tag and the mass variable. This 
weighted probability depends on the mass hypothesis and the discriminant is calculated for each 
test mass. The b-tag, the mass \ 2 an d the final discriminant distributions for all the events in 
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the low purity and high purity samples are shown in Figure for a Higgs mass hypothesis of 
105 GeV. A good agreement between the data and the expected background is observed in all 
the distributions. 

5 The HZ — >hhvi> selection 

This selection searches for events with two acoplanar jets containing b hadrons, with large 
missing energy and with missing mass consistent with ra?. 

In the first step, high multiplicity hadronic events are selected. The events are forced into 
two jets using the DURHAM algorithm. The dijet invariant mass must exceed 40 GeV. These 
requirements reduce contributions from purely leptonic two fermion final states, as well as two- 
photon interactions, while retaining a significant fraction of hadronic events from e + e~ — > qq(7) 
and W-pair production. These backgrounds are further reduced by requiring the visible mass 
to be less than 120 GeV and the mass recoiling against the hadronic system to lie between 
50 GeV and 130 GeV. 

Events from e + e~ — > qq(7) are further suppressed by requiring the longitudinal missing 
energy to be less than 0.7-^/s, the missing energy transverse to the beam axis to be greater than 

5 GeV and the missing momentum vector to be at least 16° away from this axis. The energy in 
the forward luminosity calorimeter is required to be smaller than 15 GeV. The opening angle 
between the two jets has to be greater than 69° and the angle between the dijet plane and the 
beam-axis must be greater than 3°. The b-tag distribution, after the above mentioned cuts for 
all centre of mass energies, is shown in Figure EK- 

An additional cut is applied requiring the event b-tag to be larger than 0.5. After this set 
of cuts, there are 172 events in the data, while 149 are expected from background processes. 
The efficiency for WL^hhvv with mn = 105 GeV is 62%. 

A kinematic fit imposing 4-momentum conservation and requiring the missing mass to be 
mz is performed to compute the hadronic mass. The distribution of this variable is shown 
in Figure Eb- The output of a mass independent neural network [I] is then combined with 
the hadronic mass to build the final discriminant. The distribution of the neural network 
output is shown in Figure |2fc. The final discriminant is presented in Figure 1211 for the mass 
hypothesis wh = 105 GeV. The observed data in the HZ^bbvv analysis are compatible with 
the background expectations. 

6 The HZ— >>bbe + e _ and HZ— ►bb/x + /z - selections 

The signatures for the HZ — >bbe + e~ and HZ— >bbyU + /i~ processes are a pair of high energy 
electrons or muons, with an invariant mass compatible with mz and two hadronic jets with b 
quark content. 

A high multiplicity selection is applied, also requiring two well identified electrons or muons. 
The visible energy must be larger than 0.7-^/s for the electron analysis and 0.4^ for the 
muon analysis. In the HZ— >bbe + e~ analysis, the lepton pair must have an opening angle 
greater than 100°, reduced to 90° in the HZ^bb/i + /i~ case. Moreover, in the HZ— >bbe + e~ 
channel the opening angle between the two jets must be at least 50°. The ratio between the 
transverse missing momentum and the visible energy should be less than 0.2 in the HZ — >bbe + e~ 
channel and less than 0.4 in the HZ— >bb/i + ^~ channel. The value of Y34 must be larger than 
0.0009. Finally, the invariant mass of the leptons after a kinematic fit imposing 4-momentum 
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conservation must be between 60 GeV and 110 GeV for the electrons and 50 GeV and 125 GeV 
for the muons. 

In the electron channel, 22 events are selected with 20.2 expected from the background and 
with a signal efficiency of 76% for a Higgs signal of 105 GeV. In the muon channel there are 13 
events with 9.2 expected from the background, with a signal efficiency of 56%. 

A kinematic fit that requires 4-momentum conservation and constrains the mass of the 
leptons to mi is performed. The dijet mass after this fit is combined with the b-tag values of 
the two jets, to form the final discriminant [3]. The distributions of the discriminant for the 
electron and muon channels are shown in Figure and Eh respectively, for the data compared 
to the expected background and a 105 GeV Higgs signal. The data are consistent with the 
background predictions. 

7 The HZ — >bbr + r _ and HZ — ^r + r _ qq selections 

The HZ-^bbr + r~ and HZ — >T + r~qq processes result in similar final states, partially distin- 
guished by mass and b-tag information. The semileptonic W and Z pair decays constitute the 
most significant background sources. 

Two inclusive selections are performed, one based on tau identification, the particle-based 
selection, and the other relying more on the event kinematics, the jet-based selection. Events 
are accepted if they pass either of the two selections. 

First, a high multiplicity selection is applied, also requiring a visible energy greater than 
0.4a/s. The events are required to have a value of Y34 larger than 0.0025. Background from 
e + e~ — > qq(7) is reduced by rejecting events containing photons with energies greater than 
40 GeV. The contribution of W + W _ — > qqiis (£ = e, u) is reduced by requiring the energy of 
electrons and muons to be below 40 GeV. 

In the particle-based selection, tau leptons are identified either by their decay into electrons 
or muons, or as an isolated low-multiplicity jet with 1 or 3 tracks and unit charge. In the 
jet-based selection, the event is forced into four jets using the DURHAM algorithm. Two of 
the jets must have less than 4 tracks each. These jets are considered as tau candidates, but 
at least one of them must coincide within a 3° cone with a tau candidate identified by the 
particle-based selection. Both tau candidates must be separated from the hadronic jets by at 
least 25°. Background contamination from fully hadronic W pair decays is reduced by rejecting 
events where both tau candidates decay into 3 charged particles and by requiring the visible 
energy to be smaller than 0.95^ for the particle-based selection and smaller than Q.Qyfs for 
the jet-based one. Moreover, in the jet-based selection, the missing momentum vector should 
be at least 18° from the beam axis, in order to reduce the qq(7) contamination. 

The invariant masses of the ditau and the dijet are obtained from a kinematic fit which 
imposes 4-momentum conservation. An event qualifies for the HZ^bbr + r~ channel if the 
ditau mass is consistent with the mass of the Z boson by lying between 70 GeV and 125 GeV. 
Similarly, an event is assigned to the HZ^r + T~qq channel if the dijet mass fulfils the same 
requirement. The opening angle of the taus or jets assigned to the Higgs boson must be 
larger than 70°. Those assigned to the Z must be at least 100° apart. Events selected by the 
HZ— >bbe + e~ and the HZ— >bbu. + u,~ analyses amount to 3% and are not considered here. 

In total, 19 events pass either of the tau selections, with 22.3 events expected from back- 
ground processes. The efficiency is 29% for both HZ— >bbr + r~ and HZ— »T + T~qq at mn = 
105 GeV. 
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The discriminant is defined as in the HZ — >bbe + e~ and HZ—>-hh[i + {i~ analyses. Events 
passing both HZ — >bbr + r~ and HZ^T + r~qq selections are assigned to the channel for which 
the value of this discriminant is higher. This discriminant is used as the final variable for the 
HZ^bbr + T~ selection. For the HZ-^T + r~qq selection, the mass of the tau pair, calculated 
by constraining the invariant mass of the two other jets to mz, is used as the final discriminant. 
Distributions of the final discriminants are shown in the Figures Et and Eli for data, expected 
background and a 105 GeV Higgs signal. No evidence of a signal appears in either of the tau 
channels. 

8 Combined results 

The results of all the previously described analyses are combined together to set a lower limit 
on the mass of the Standard Model Higgs boson. A combined CL on the absence of a signal is 
derived from the distributions of final discriminants in a scan over ma from 50 GeV to 110 GeV. 
The CL is calculated using the technique of References [H3HI], which takes into account the 
correlated and the statistical errors. 

The systematic uncertainties on the signal and background expectations are derived using 
the same procedure adopted in previous Standard Model Higgs searches [TH] . The overall sys- 
tematic uncertainty is estimated to be 10% on the number of background events and 4% on 
the number of signal events. The statistical uncertainty on the background arising from the 
finite number of generated Monte Carlo events is uncorrelated from bin to bin in the final dis- 
criminant distributions, and has little effect on the CL. Bins with a signal-to-background ratio 
below 0.05 are not considered in the calculation of the CL. This cut was chosen to maximise the 
median CL, as calculated from a large number of Monte Carlo experiments, thereby minimising 
the degradation of the result due to these systematic and statistical uncertainties. The results 
of all the analyses after such a signal-to-background cut are summarised in Table ^for the data, 
the expected background and for Higgs signals of 100 GeV and 105 GeV. The number of signal 
events includes cross-efficiencies from other channels, fusion processes and charm and gluonic 
Higgs decays. 

The measured value of the CL as a function of the Standard Model Higgs boson mass in the 
range 95 < mn < 110 GeV is shown in the Figure HJi, along with the median of the CL distribu- 
tion as calculated from a large sample of Monte Carlo experiments under the background-only 
hypothesis. The median CL represents the sensitivity of the analysis and is equal to 95% at 
m H = 105.2 GeV. The results of L3 Standard Model Higgs searches at lower centre of mass 
energies j3] are included in the calculation of these confidence levels. Values of run from 50 GeV 
to 95 GeV are excluded in the Standard Model with a confidence level greater than 99%. 

Figure Eb shows the background confidence level, CLf,. The observed CLj, is the probability 
to observe a smaller number of events than the one actually observed, for a background-only 
hypothesis. The expected CL& for a background-only hypothesis is 0.5. The observed CL& gives 
a measure of the consistency of the data with the expected background. The relatively high 
value of CL;, in the mass region between 96 and 100 GeV reflects a 2.2 sigma excess of data 
relative to the background predictions in this mass region. This excess is mainly due to few 
candidates in the HZ— >bbqq channel and to one candidate in the HZ— >bbqq channel with a 
reconstructed Higgs mass of 100.5 GeV. For mn =107 GeV, where the observed CL falls below 
95%, the CL b is 8%. 
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The lower limit on the Standard Model Higgs boson mass is set at 

m H > 107.0 GeV at 95% CL. 

This new lower limit improves upon and supersedes our previously published results. A similar 
result from the 192-202 GeV data was reported ^H] • Results from the year 2000 LEP run were 
recently published [THIGH! • 
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Table 1: The number of observed candidates (N D ), expected background events (N B ) and 
expected signal (N s ) for the = 192 GeV, 196 GeV, 200 GeV, 202 GeV data after a cut on 
the final discriminant corresponding to a signal-to-background ratio greater than 0.05. This 
cut is used to calculate the confidence levels. 
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Figure 1: The b-tag (a and b), the mass % 2 (c and d) and the final discriminant (e and f) 
distributions for the low purity and high purity HZ — >bbqq selections. The points are the 192- 
202 GeV data, the open histograms are the expected background and the hatched histograms 
are the expected Higgs signal from the HZ — >bbqq channel with mu = 105 GeV. The last bin 
on the right of the b-tag distributions groups the overflows. 
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Figure 2: Distributions of (a) b-tag, (b) hadronic mass (c) neural network output and (d) final 
discriminant for the HZ — >Yfovv selection. The points are the 192-202 GeV data, the open 
histograms are the expected background and the hatched histograms the HZ — >bbvi? expected 
signal with mu = 105 GeV. 
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Figure 3: Distributions of the final discriminant for the (a) HZ^bbe + e~, (b) HZ^bb/z + /i~, 
(c) HZ^bbr + T~ and (d) HZ^r + r _ qq selections. The 192-202 GeV data (points) are com- 
pared to the expected background (open histogram) and to the expected Higgs signal (hatched 
histogram) of 105 GeV. The signal events in the HZ— »bbr + T~ and HZ-^T + r~qq distributions 
include the cross-efficiencies for these channels. Events are uniquely assigned to only one of 
these channels. 
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Figure 4: (a) The observed (solid line) and expected median (dashed line) signal confidence 
levels (1-CL) as a function of the Higgs mass, (b) the observed (solid line) and expected (dashed 
line) background confidence level CL5 as a function of the Higgs mass. The lower limit on the 
Higgs mass is set at 107.0 GeV at the 95% CL. 
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